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ABSTRACT 



Context. Relativistic effects in the oscillatory damping of magnetic disturbances near two-dimensional X-points are 
investigated. 

Aims. By taking into account displacement current, we study new features of extremely magnetized systems, in which 
the Alfven velocity is almost the speed of light. 

Methods. The frequencies of the least-damped mode are calculated using linearized relativistic MHD equations for wide 
ranges of the Lundquist number S and the magnetization parameter a. 

Results. The oscillation and decay times depend logarithmically on S in the low resistive limit. This logarithmic scaling 
is the same as that for nonrelativistic dynamics, but the coefficient becomes small as ~ <t~ 1//2 with increasing a. These 
timescales approach constant values in the large resistive limit: the oscillation time becomes a few times the light 
crossing time, irrespective of a, and the decay time is proportional to a and therefore is longer for a highly magnetized 
system. 
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1. Introduction 

The importance of magnetic reconnection manifests itself in 
various energetic astrophysical phenomena, including rela- 
tivistic objects such as pulsars, magnetars, active galactic 
nuclei and gamma ray bursts. The characteristic propaga- 
tion velocity for magnetic disturbances, the Alfven velocity, 
depends on the magnetization parameter a: 2a represents 
the ratio of the magnetic to the rest mass energy density 
of the plasma. When the magnetization parameter (j>l, 
the Alfven velocity is almost the speed of light; the veloc- 
ity becomes nonrelativistic in the opposite limit, a <C 1. In 
this paper, we consider some inherently relativistic features 
that may appear in the magnetic reconnection when a is 
large. 

In a simple analysis of the Sweet-Parker type recon- 
nection, the structure of the reconnection layer depends 
on two large dimensionless numbers: a and the Lundquist 
(or magnetic Reynolds) number S, an inverse of resistiv- 
ity (|Lvutikov fc Uzdenskvl (|2003h ). For a < S, the in- 
flow velocity is nonrelativistic, and the reconnection is very 
similar to the classical Sweet-Parker model. However, for 
a ^> S 1, the inflow velocity becomes relativistic. 
iLvubarsk v (2005) incorporated the compressibility of mat- 
ter and found that the inflow velocity is always sub-Alfvenic 
an d remains much less than the speed of light, contradict- 
ing iLvutikov fe Uzdenskvl (|2003l) . The reconnection rate is 
still estimated by substituting c for the Alfven velocity in 
the nonrelativistic formula, even in the relativistic regime. 

Small differences may also originate from the as- 
sumption of a steady state. Numerical simulations 
of an anti-parallel magnetic configuration in two di- 
mensions have been performed without assuming a 
steady state using re l ativist ic resistive MHD code 
(jWatanabe fc: Yokovamal (|2006l) ). a relativistic two- fluid 



model lfZenitani et al. (120091)). an d PIC simulations on ki- 
neti c scale (IZenitani fc Hoshind (|2005l l2007t 12008(1 ). Sec 
also iKomissarovl (|2007|) for the numerical schemes of the 
resistive relativistic MHD. These approaches have clearly 
demonstrated the relativistic dynamics, but simulation in 
a wide range of parameters would be time-consuming. 
Moreover, the resolution becomes poor for small resistiv- 
ity. 

The dynamics at an X-type null point, where a cur- 
rent sheet forms and the magnetic energy is dissipated, 
have been stu died previously. In c o ntext of nonrelativis- 
tic dynamics, ICraig fc McClvmontl (|1991l) considered the 
behavior of MHD waves near the X-point in the cold 
plasma approximation using linear perturbation theory. 
They showed the remarkable result that the dissipation 
time behaves as ~ (ln^ 1 ) 2 . This logarithmic dependence, 
in contrast to the normal power behavior ~ S a , indi- 
cates fast d ecay. Subseque ntly, the problem was studied 
analvticallv (|Hassaml (|1992D ) and by considering the prq p- 
agation of linearized waves dMcLaughlin fc Hoodl (l200l ). 
Some physical properties of a more realistic system have 
also bee n included, such as no n - linear waves with thermal 
pressu re (jMcClvmont fc Craid (I1996D: iMcLauehlin et all 
(|2009D ). elect ron inertial effects dMcClements et al.l (|2004l )). 
and viscositv (|Craig et al.l (I2005D: ICraid (I2008D). Se e a re- 
cent review of this topic given bv lMcLauehlin et al.l (|2010l) 
and references therein. 

The main concern of this paper is to explore relativis- 
tic effects on the dynamical reconnection at an X-point. In 
particular we consider whether the reconnection is qualita- 
tively modified for a highly magnetized system with a ^> 1, 
such as a magnetar. We adopt a very simple system in or- 
der to understand the diff erences, if any. Our wor k is a rel- 
ativistic extension of ICraig fc McClvmon t (1991). That is, 
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we will calculate complex normal frequencies, which deter- 
mine the oscillatory damping of the magnetic disturbances 
with small amplitudes, neglecting thermal pressure, viscos- 
ity and so on. The problem may be solved as an initial value 
problem, but the initial data inevitably contain electromag- 
netic waves besides MHD waves, and subsequent evolution 
may be complex. In section 2, we discuss our numerical 
methods and boundary conditions. Our results are shown 
in section 3. Section 4 contains our conclusions. 



2. Model 

2.1. Equations for relativistic dynamics 

We consider a two-dimensional problem, assuming d/Bz — 
0. In our model, the magnetic held B is located on a plane 
and the electric held is perpendicular to it, E = Ee z . The 
electric current je z is also perpendicular to the plane, and 
the charge density consistently vanishes, since V E = 0. 
These electromagnetic fields can be expressed in terms of 
only the z-component of a vector potential A — Ae z as 



VA x e 2 , 



E = — 



IdA 
c dt' 



(1) 



The flux function A satisfies with a wave equation with a 
source term: 



c 2 dt 2+V I- 1 - 



(2) 



where the displacement current is included in contrast to 
the usual nonrelativistic treatment. 

The dynamics of the plasma flow is determined by the 
continuity equation 



dp 



+ V • (pv) = 0, 



and the momentum equation with the Lorentz force 

p + v ■ V ) jv = ] -e z x B = 3 -VA, 
\ot c c 



(3) 



(4) 



where 7 = (1 — (u/c) 2 ) -1 / 2 , p is the mass number den- 
sity in the laboratory frame, and the proper one is p/j. 
In eq. ([!]), the Coulomb force vanishes and thermal effects 
in the pressure and internal energy are neglected in the 
cold limit. This cold plasma approximation simplifies the 
problem:The slow magnetoacoustic wave is absent. In non- 
relativistic dynamics, it is found that propagation of the 
fast one causes the current density to accumulate at the X 
point , where the energy i s dissip ated (|McLaughlin fe Hoodl 
12004 iMcLaughlin et all . l2009tl . Thermal pressure is ne- 
glected, since our concern is the propagation in linearized 
system. The finite pressure is meaningful in fully non-linear 
dynamics, where coupling and mode conversion between 
MHD waves are important in the neighborhood of the dis- 
sipation zone. 

Ohm's law with resistivity 77 can be written as 



E 



1 



(v x B). 

C JC 



—3> 



(5) 



The relativistic motion redu ces the resistivity by th e 
Lorentz factor 7. (See, e. g., iBlackman fc Field! (|l993h : 
iLvutikov fc Uzdenskvl (|2003f) .i However, this factor may be 
set to 7 = 1 for a linear perturbation from a static back- 
ground. 

2.2. Normal mode for the linearized system 

We consider the dynamics of small perturbation in the 
vicinity of an X-point, which is governed by current-free 
(Jo = 0), static (vq = 0) background fields with uniform 
density (p = p ). 

The magnetic potential Aq of the background field can 
be written in the Cartesian (x, y) or polar coordinates (r, 9) 
as 



(7) 



where L is a normalization constant for the length and Bq 
is a constant representing the magnetic field at r = L. 

The linear perturbation approximation for eqs. (H])-® 
reduces to a single equation for 5 A: 



d , (VA ) 2 



dt 47rpo 



1 f) 2 \ B 2 

-#d¥ + v ) SA -W> 6A = - (8) 



By using normalized length r = r/L and time t = vot/L, 
where vq — Bq/IAttpq) 1 ^ 2 , eq. §E§ becomes 



1 d 



B 2 - \ B 2 
<r^+V 2 )5A-^5A = 0, (9) 



^s, dt ' ' J \ dt 2 ) dt 2 
where s* and a are non-dimensional parameters given by 



(10) 



Attpqc 2 



(11) 



The magnetization parameter a has been introduced 
through the displa cement current, and h ence eq. (jSJ) be- 
comes eq. (2.4) of ICraig fc McClvmont] (fl99l when the 

D'Alembertian — cjp- + V 2 is replaced by the Laplacian V 2 
in the limit of a = 0. It should be noted that vq represents 
the Alfven velocity at radius L only in the nonrelativistic 
case. The Alfven velocity at L is in general given by Va = 
ca 1/2 /(<r + 1) 1/2 = v /(<7 + 1) 1/2 . For highly magnetized 
cases where a 3> 1, we have Va ~ c, whereas Va ~ «o for 
J < 1. Although eq. © is used for mathematical calcula- 
tion, the physical results are presented after normalization 
by Va- The Lundquist number S characterizing the system 
is defined in terms of the Alfven velocity Va, the radius L 
and resistivity r\ as 



S = 



V A L 



(12) 



which, in terms of A, is 



The related parameter s* is s* = (cr + 1) 1 / 2 S'. 

Equation (O exhibits two different behaviors near to 
and far from the origin. For large r, the dissipating term 
with s» can be neglected, so that we have 
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This is exactly the equation in the cold plasma limit for the 
propagation of a fast magnetoacoustic wave, whose velocity 
at f is given by the Alfven velocity 



v A (r) 



v Q r 



CO 



1/2; 



(af 2 + l)i/ 2 (af 2 + l)i/ 2 ' 



(14) 



On the other hand, close to the origin, the term with f 2 
can be neglected in eq. ©. After integrating by t once, we 
have 



'dp 



l_d_ 

s» di 



+ V 



SA = 0. 



(15) 



This is the so-called telegraphist's equation, in which the ef- 
fect of the finiteness of the velocity c on the resistive losses, 
or the effect of resistiv ity on the wave eq u ation , is taken 
into account. (See, e.g. lMorse fc Feshbacb) (fl953lU In the 
limit of a = 0, the equation becomes the diffusion equa- 
tion. Thus, eq. (|9]) leads to an advection-dominated outer 
region described by eq. (|13p and a diffusion dominated in- 
ner one described by eq. dTB]) . The diffusion region may be 
highly modified in nature for large c, as electromagnetic 
wave propagation becomes important even in the diffusion 
zone for a highly magnetized system. The critical radius f c , 
which separates the two regions, will be determined by the 
following normal mode analysis. 

We solve eq. (O as an eigenvalue problem in the form 

SA = /(f) cxp(im0) exp(-iCot) 

= f(f)exp(im9)exp(-iu)VA(a + l) 1/2 t/L), (16) 

where w is a complex number. We only consider the axially 
symmetric m = mode, whi ch is relevant to reconnectio n 
at the origin, as discussed in ICraig fc McClvmontl (|1991l) . 
Another type of re conn ection for m ^ is dis cussed by 
lOfman et all (|l993l) and lVekstein fc Bianl (|2005h . but that 
not is considered here. Equation © becomes 



1 d _ d , 
——r — j 

r dr dr 



f = 0. 



(17) 



From this, a natural choice of the core radius r c is of or- 
der (\Co\l s-t) 1 / 2 ~ iS*" 1 / 2 and f c c orresp onds to the usual 
skin depth ()Craig fc McClvmontl (|l99lh ). The dissipative 
term is dominant for f < f c , whereas outside the critical 
radius eq. ([!?)) represents wave propagation, since the term 



with 



can be neglected. The current density is 



concentrated around the null point. 

A series solution inside the radius f c may be expressed 

as 



/ = 1 - ~ {Co 2 o + is*Co)f 2 + 



(18) 



where we have normalized to / = 1 at the origin. We solve 
eq. (|17[) with boundary condition (|18[) . from f = f c to 1, 
assuming a complex number Co. The boundary condition 
imposed on the circle f = 1 is / = 0. This means that the 
magnetic flux is frozen and SE = Sj = 5v = there. Thus, 
we have a one-dimensional eigenvalue problem for Co. 

Our main concern is not whole eigenfrequency spec- 
trum, but rather the lowest frequency mode, which persists 
for a long time in the magnetic reconnection. In particu- 
lar, we will study the effect of the magnetization parameter 
on it. For this purpose, we first calculate Co for the case 
cr = 0, and then repeat the calculation, gradually changing 
the parameter S or a . 



3. Results 

The oscillation time t osc is defined in terms of the real part 
of the eigenfrequency Co by t osc — 2nL/ (((T+l) 1 / 2 Re((D) Va). 
(A factor (er+1) 1 / 2 comes from our normalization of Co. (See 
eq. (IT6l) .) Figure 1 shows the normali zed time VAt nsr JL as 
a func tion of S for several values of a. ICraig fc McClvmontl 
(|!99lh showed that the relation V A t sc/L « 2 In 5 « 
4.6 log S holds for a wide range of S with a — 0. The ori- 
gin of this relation can be understood by considering the 
traveling time of an MHD wave from the outer boundary 
to the resistive region, 



L 



-dr. 



r. VA(r 



(19) 



The velocity in the limit of a = is scaled by va oc f , and 
the dominant contribution in eq. (|19[) comes from a small 
core region. By choosing the lower boundary f* as f c , we 



have t 



lnr c oc In S. 



When a is included, the oscillation time deviates from 
the relation VAt osc /L ~ 21ns*. The normalized time, in 
general, becomes smaller than that at a — 0, as shown 
in Fig. 1. The logarithmic dependence with S can be seen 
only in the larger regime, and the coefficient in front of 
In S becomes smaller as a increases. The Alfven veloc- 
ity becomes relativistic for a > 1 at the boundary, and 
approaches zero toward the center. The velocity becomes 
nonrelativistic, at the radius fjv ~ er -1 / 2 for a 3> 1, and 
the velocity is almost equal to c outside this radius. The 
wave traveling time in eq. (I19[) is almost determined by the 
slow region inside f jy , and the system size may be regarded 
as being effectively reduced to ct _1 / 2 L. We therefore have 
VkW(^~ 1/2 £) w 21n5, i.e, V A t osc /L « 2a- 1 / 2 In S for 
the large S regime. This property can be seen from the 
curves around log S ~ 50 in Fig. 1, except for a = 10 4 . A 
factor of (a + l) 1 / 2 instead of a 1 / 2 may provide a better 
extension to a = 0, but a simple correction is used here. 

Figure 1 also shows that VAt osc /L approaches a con- 
stant in the small S regime, for sufficiently large a. 
Asymptotically the value of this constant as S — > 1 is em- 
pirically V A t osc /L ~ ct osc /L w 2.5, which is independent of 
a, as far as a > 10 2 . In our model, the core size increases 
as r c oc S" -1 / 2 , and hence the traveling time (fT9")) becomes 
smaller with decreasing 5, but the lower bound is a few 
times the light crossing time for a region of size L. 

The critical value S c , which discriminates between con- 
stant VAtosc/L for smaller S and VAt OS c/L oc lnS* for larger 
S, is given approximately by lnS* c ~ cr 1 / 2 , or logS'c ~ 
0.4CT 1 / 2 . The transition is not very sharp but the relation 
does give the approximate boundary between two distinct 
behaviors. Because logS'c ^ 1 for cr = 10 and logS'c ~ 40 
for a = 10 4 , which are located at the edges of Fig. 1, the two 
different behaviors are not clearly shown for these param- 
eters. This critical value S c also characterizes a transition 
in the decay time as will be discussed below. 

The decay time is related to the imaginary part of Co, 
tdecay = L / {{o + 1) x l 2 |Im(w) | V A ) ■ Figure 2 shows the nor- 
malized decay time V At decay I 'L as a function of S for sev- 
eral values of a . The time for a = scales as VAtdecay / L 
= 2(lnS) 2 /7r 2 (ICraig fc McClvmontl (119911^ . This scaling 
relation is also broken by the inclusion of a. The small 
and large S regimes are different, as they are for the os- 
cillation time. A typical example is given by the curve for 
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Fig. 1. Normalized oscillation time t osc = VAt OS c/L as a 
function of Lundquist number 5, for magnetization param- 
eter values a = 0, 10, 10 2 , 10 3 , 10 3 5 and 10 4 . 
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Fig. 2. Normalized decay time r decay = VAtdecay/L as a 
function of Lundquist number S 1 , for magnetization param- 
eter values a = 0, 10, 10 2 , 10 3 , 10 35 and 10 4 . 



a = 10 3 : the critical value is log S c ~ OAa 1 / 2 ~ 13 for this 
case. Logarithmic dependence can be seen for logS 1 > 20, 
whereas the curve becomes constant for logS" < 7. The 
relation VAtdecay/L oc (lnS*) 2 can be seen in the large S 
regime, S 3> S C) except for a = 10 4 , but the timescale is 
reduced to approximately VAtdecay/L s» 2er~ 1 / 2 (ln S) 2 /n 2 
for a 3> 1. The factor c -1 / 2 can be interpreted as being 
due to an effective reduction of the system's size, as con- 
sidered for the oscillation time. The normalized decay time 
becomes the minimum around S c . 

In the small S regime, S <C S c , normalized decay time 
approaches a constant value VAtdecay/L ~ 0.14c. The nor- 
malized decay time for fixed S increases with the magneti- 
zation parameter a. The limit of a — > oo corresponds to the 
vacuum, in which there is no matter (p = 0) and the dis- 
sipation time becomes infinite. This cr-dependence comes 
from taking account of the finiteness of c in the resistive 
losses. (See eq. (fT5)) .) This effect can be neglected in the 
large S regime, where the approximation of instantaneous 
dissipation is good. However, the effect becomes evident in 
the small S regime. 

The energy E of perturbation decreases due to the 
Ohmic dissipation 



dE 
~dt 



-V / fdV. 



(20) 



The linearized form with Fourier component provides an 
expression of the decay time as 

y^decay _ ^Jq {^B + Se E + Se M^TTfdf 



Jo \Sj\ 2 2Trrdr 



(21) 



where 5s is dimensionless energy density of magnetic field, 
electric field, kinetic energy of the fluid, and Sj is dimen- 
sionless current density. Their explicit forms are given by 



Set 



1 



\5E\- 



a 



5e 



M 



:Po\Sv\ 



2\Q 



r r ar ar 



and 



(22) 
(23) 
(24) 

(25) 



Spatial distributions of these energy densities are displayed 
in Fig. 3 for S = 10 5 , a = 10 1 and in Fig. 4 for S = 10 5 , 
(7 = 10 4 . These functions are calculated by numerical so- 
lution outside f c , and by the analytic asymptotic form 
eq. (fT5|) inside it. Note that a sharp peak in Sem and 
8e~B is located within f c . Both kinetic energy of matter 
and magnetic energy are accumulated from outer part to 
the core(~ f c ), and are dissipated in the central region. 
However, distribution of electric energy is flat. These overall 
features are not so much different in Figs. 3 and 4, although 
the sharp peak shifts by f c = (w/((a + l) 1 / 2 ^)) 1 / 2 . 

The magnitude of See is much smaller than that of Seb 
in Fig. 3 (a = 10 1 ), whereas See becomes comparable to 
Seb in Fig. 4 (cr = 10 4 ). The electric energy is approxi- 
mately proportional to cr, as shown in ea. (|2"3"|) . and signif- 
icantly contributes to the sum of energy. Hence, the de- 
cay time becomes longer with the increase of cr for fixed 
S, since the total energy increases. (See ea. (|2"Tj) .) In the 
large S regime, however, the functions Seb and Sem are 
much larger than See, so that the electric energy can be 
neglected. The decay time does not increase with a in this 
regime. 

4. Discussion and conclusions 

Relativistic MHD differs, in general, from the nonrelativis- 
tic case in at least three ways: (i) the Lorentz factor 7, 
(ii) the Coulomb force p e E, and (iii) the displacement cur- 
rent c~ l dE/dt in Maxwell's equation. The Lorentz factor 
appears in the flow velocity and also in the resistivity of 
Ohm's law as a Lorentz contraction. The difference is of 
order (w/c) 2 in magnitude. Since we considered a linear 
perturbation from the static state, the inflow velocity is 
not very large and the Lorentz factor may approximate to 
7 = 1. The magnitude of p e E is of order (v/c) 2 times the 
Lorentz force j x B, and is hence neglected in nonrelativis- 
tic MHD. Moreover, the charge density is always zero due 
to the 2D X-point geometry considered here, so that the 
Coulomb force p e E vanishes exactly. This leaves the dis- 
placement current as a possible factor for the difference be- 
tween relativistic and nonrelativistic MHD. We have stud- 
ied its effects, especially on the dynamics of the magnetic 
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Fig. 3. Normalized energy density Se as a function of x = 
lnf for S — 10 5 and a — 10 1 . The function Sem has a sharp 
peak, and is shown with a reduction factor 4 x 10 -25 , while 
See is magnified by 2 x 10 3 . 
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Fig. 4. Normalized energy density Se as a function of x = 
lnf for S = 10 5 and c = 10 4 . The function Sem is shown 
with a factor 6 x 10~ 14 , while See is shown with a factor 6. 



reconnection using a simplified system based on linearized 
equations in the cold plasma limit. The magnetization pa- 
rameter a is incorporated in the basic equation through the 
displacement current and the oscillation and decay times 
for the least-damped mode were calculated numerically for 
parameters S =10-10 50 and a = 0-10 4 . 

In the system with a = 0, for which the displacement 
current can be neglected, the oscillation and decay times are 
proportional to In S and (lnS 1 ) 2 , respectively. By including 
a, these timescales are modified in different ways, in two 
regimes, which are characterized by S ^ S c or S <C S c for 
S c ?» exp(cr 1 / 2 ). For low resistivity, S 3> S c , a logarithmic 
dependence with S can seen, but the timescales normalized 
by the boundary radius L and the Alfven velocity Va be- 
come smaller with increasing a. The smaller timescales can 
be explained as being due to an effective reduction in the 
size of the system, or the enlargement of the outer region 
where MHD waves propagate at almost the speed of light 
and the traveling time is negligible. On the other hand, for 
high resistivity, S <C S c , a new feature appears in both 
the oscillation and decay times, which do not depend on S. 



The oscillation time is a few times the light crossing time 
and does not depend on ex. The dissipation time becomes 
longer in proportion to a and goes to infinity in the limit of 
<7 — > oo, that is, no dissipation in the vacuum. Reconnection 
at the X point is thought to be "fast" , since the dissipation 
time is scaled with (lnS 1 ) 2 . Actual time is of the order of 
10-10 3 times crossing time with Alfven velocity. The dis- 
placement current significantly spoils the good property, 
and the timescale increases with a in high resistive region. 
The increase of the decay time is related with deficiency of 
matter, which is involved in the Ohmic dissipation. 

Magnetic reconnection is expected to be an important 
process of abrupt energy release in the solar and magne- 
tar flares. For example, the explosive tearing-mode recon- 
nection in the magnetar like the sola r flares is discussed 
(jLvutikovl (|2006l) ; lMasada et all (|2010D ). Dimensionless pa- 
rameters are however quite different in them: a ~ 10~ 4 
and S ~ 10 14 in solar corona, whereas it is likely that 
(T>1 and S 3> 1 in a magnetar magnetosphere. Present 
result in an X-type collapse suggests the dissipation time 
t ~ O.IctL/Va ~ 10~ 5 cr (L/10 6 cm) s under highly magne- 
tized environment. The spiky rise time (< 0.1s) or short 
duration (< Is) of the magnetar flare may significantly 
constrain ah. The energy of the flare AE(~ 10 45 erg) 
should be a part of magnetic energy within the volume L 3 : 
BqL 3 ~ po&L 3 > AE. These two conditions provide an 
upper limit of a as a < 10 4 5 (po/(g/cm )) x ^ 2 . In such high 
energy events, radiation and possibly pair creation may be 
important in the energy transfer. Further study is needed 
for these effects. However, the results in this paper demon- 
strate that the dynamics significantly depends on the mag- 
netization parameter through the displacement current. 
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